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Abstract. I present results from numerical simulations of gas dynamics outside luminous accre- 
tion disks in active galactic nuclei. The gas, gravitationally captured by a super massive black 
hole, can be driven away by the energy and momentum of the radiation emitted during black 
hole accretion. Assuming axisymmetry, I study how the mass accretion and outflow rates, and 
the flow dynamics respond to changes in radiation heating relative to radiation pressure. I find 
that for a 10 8 Mq black hole with the accretion luminosity of 0.6 of the Eddington luminosity 
the flow settles into a steady state and has two components: (1) an equatorial inflow and (2) a 
bipolar inflow/outflow with the outflow leaving the system along the disk rotational axis. The 
inflow is a realization of a Bondi-like accretion flow. The second component is an example of a 
non-radial accretion flow which becomes an outflow once it is pushed close to the rotational axis 
where thermal expansion and radiation pressure accelerate it outward. The main result of this 
preliminary work is that although the above two-component solution is robust, its properties 
are sensitive to the geometry and spectral energy distribution of the radiation field. 
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1. Introduction 

The radiation properties of active galactic nuclei (AGN) and the AGN central location 
in their host galaxies imply that they play a very important role in determining the 
ionization structure and dynamics of matter not only in their vicinity but also on larger, 
galactic and even intergalactic scales (Cotti & Ostriker, 1997, 2001; King 2003; Murray, 
Quataert, & Thompson 2005; Sazonov et al. 2005; Springel, Di Matteo & Hernquist 
2005; Hopkins et al. 2005, and references therin). Many observational results support 
this suggestion, in particular the presence of broad emission and absorption lines in 
AGN spectra. The ionization structure and dynamics of the gas responsible for these 
lines can be driven by radiation, even for sub-Eddington sources. The driving can be due 
to radiation pressure or radiation heating, or both (e.g., Begelman, McKee and Shields, 
1982; Shlosman, Vitello & Shaviv 1985; Ostriker, McKee, & Klein 1991; Arav & Li 1994; 
Murray et al. 1995; Proga, Stone & Kallman 2000; Proga & Kallman 2002, 2004). 

In this paper, I present results from hydrodynamical simulations of a non-rotating gas 
on sub-parsec- and parsec-scales in AGNs. I use a simplified version of the numerical 
method developed by Proga, Stone & Kallman (2000) to study a related problem of ra- 
diation driven disk winds in AGN (for details see Proga et al. 2006 in preparation). I 
consider an axisymmetric flow accreting onto a supermassive black hole (BH). The flow 
is non-spherical because it is irradiated by an accretion disk. The disk radiation flux is 
the highest along the disk rotational axis and is gradually decreasing with increasing 
polar angle, 6 as cos#. The flow is also irradiated by an isotropic corona (see eq. 2.1). I 
take into account the radiation heating and cooling, radiation pressure due to the elec- 
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tron scattering and spectral lines. I adopt a simplified treatment of photoionization, and 
radiative cooling and heating allowing for a self-consistent calculation of the ionization 
state, and therefore the line force, in the flow. 

2. Results 

I assume the mass of the non-rotating BH, Mbh = 10 8 M and the disk inner 
radius, r* = 3 rs = 8.8 x 10 13 cm throughout this paper. I consider the case 
with the rest mass conversion efficient r] = 0.0833 and the mass accretion rate, M a 
= 10 26 g s _1 (= 1.6 M0yr -1 ). These system parameters yield the accretion luminosity, 
L = 7.5 x 10 45 erg s _1 , corresponding to 0.6 of the Eddington luminosity. To determine 
the radiation field, I specify the fraction of L in the UV and X-ray band, as /uv and /x, 
respectively. 

Although, a quasar is powered by accretion, the disk accretion rate, that determines 
radiation, is not coupled to the rate, Mm(n) at which the non-rotating gas leaves the 
computational domain of these simulations through the inner boundary. This M a -Mi n (ri) 
decoupling is physically motivated because an accretion disk is build from rotating gas 
but non-rotating gas will not significantly contribute to the disk mass and to the system 
luminosity. In this sense, I study radiatively inefficient flows accreting onto an object 
with a radiatively efficient accretion disk. 

I present here preliminary results from simulations where all model parameters are 
fixed except for /uv and fx- I consider three cases: case A with /uv = 0.5 and fx = 0.5, 
case B with /uv = 0.8 and fx — 0.2, and case C with /uv = 0.95 and fx = 0.05 (see 
Table 1 for summary of the runs). 

At large radii, the radial radiation force from the disk and spherical corona can be 
approximated as 

(r, 0) = £^ [2 cos0/ uv (l + M(t)) + fx] , (2.1) 

where r and 6 are the radius and polar angle in the spherical polar coordinate system, 
respectively while the terms in the brackets with f\jy and fx correspond to the disk 
and corona contribution, respectively. For simplicity, I assume that all UV photons are 
emitted by the disk whereas all X-rays are emitted by the corona. The M(t) is the so- 
called force multiplier - the numerical factor which parameterizes by how much spectral 
lines increase the scattering coefficient compared to the electron scattering coefficient 
(Castor, Abbott & Klein 1975). 

Figure 1 compares the results from run A and run C. I specified the outer boundary 
in the following way. The density and temperature at the outer radius, r Q was set to 
Po = 10 -21 g cm -3 and Tq = 2 x 10 7 K, respectively. The velocity was set to zero. At 
the outer radial boundary, during the evolution of each model I continue to apply these 
constraints that the density and temperature are fixed at constant values at all times. For 
the initial conditions, I set all variables constant and equal to their values at the outer 
boundary, as listed above. The figure shows the instantaneous density, temperature, and 
distributions and the poloidal velocity field of the models. Additionally, it also shows the 
so-called Compton radius corrected for the effects of radiation pressure due to electron 
scattering 

R c =Rc[l~ r(2cos0/ uv + fx)], (2.2) 
where Rq = GM^nm p / kTc = 8.03 x 10 18 cm =9.1 x 10 4 r* is the uncorrected Compton 
radius for the Compton temperature, Tq = 2 x 10 7 K. 

In all runs, the flow settles quickly into a steady state (within ~ a few xlO 12 s which 
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Figure 1. Comparison of the results for run A and C (left and and right column, respectively). 
Top row of panels: Maps of logarithmic density overplotted by the direction of the poloidal 
velocity. Bottom row of panels: Maps of logarithmic temperature overplotted by the direction 
of the poloidal velocity. The solid curve in the bottom left corner marks the position of the 
Compton radius corrected for the effects of radiation pressure due to electron scattering (see eq. 
2.2 in the main text). The length scale is in units of the inner disk radius (i.e., r' — r/r* and 
z — z/r*). The computational domain is defined to occupy the angular range 0° ^ ^ 90° and 
the radial range n = 500 r* ^ r ^ r Q = 2.5 x 10 5 r*. 



correspond to a few dynamical time scales at r Q , r = (^/GMbh) 1 ^ 2 = 9 x 10 11 s ). 
The steady state consists of two flow components (1) an equatorial inflow and (2) a 
bipolar inflow/outflow with the outflow leaving the system along the pole. The outflow 
is collimated by the inflow. The calculations capture the subsonic and supersonic parts 
of both the inflow and outflow. Although the same components can be identified in all 
runs, their size, density and temperature, and the degree of outflow collimation depend 
on /uv and fx (the spectral energy distribution and geometry of the radiation field). 
In particular, the outflow power (e.g., measured as the kinetic energy, Pk, and thermal 
energy P t h carried by the outflow) and the degree of collimation is higher for the model 
with the radiation dominated by the UV/disk emission (run C) than for the model 
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Table 1. Summary of results 



Run 


/uv 


fx 


Min(r ) 


Min(n) 


Mout(r ) 


V r 


Pk(r ) 


Pt(To) 


A 


0.5 


0.5 


4 


1 


3 


700 


2 


4 


B 


0.8 


0.2 


8 


3 


5 


4000 


100 


1 


C 


0.95 


0.05 


9 


1 


8 


6700 


300 


0.03 



The quantities in the table are in the following units: Mi n (r ), Mi n (n), and M ou t(r ) are in 
units of 10 25 g s _1 , v r is in units of km s _1 , and Pk(r ) and Pt(r ) are in units of 

10 10 erg s" 1 . 

with the radiation dominated by the X-ray /corona emission (run A). I note that a very 
narrow outflow driven by radiation pressure on lines can carry more energy and mass 
than a broad outflow driven by thermal expansion (compared results for runs C and A). 

3. Conclusions 

The simulations show that AGN can have a substantial outflow which originates from 
the inflow at large radii. Such an outflow can control the rate at which non-rotating matter 
is supplied to the AGN central engine because the outflow mass loss rate, M ut(n), can 
be significantly higher than the mass inflow rate at small radii, Mm(n)- For example, in 
run C, as little as 10% of the inflow at large radii reaches small radii because 90% of 
the inflow is turned into an outflow. However, even the power of the strongest outflow is 
very low compared to the radiation power (i.e., for run C, P^jL — 4 x 10 -4 ). Finally, the 
inflows and outflows, found in these simulations, can be related to material responsible for 
broad absorption and emission lines, and narrow absorption and emission lines observed 
in X-ray and UV spectra of AGN. 
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